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Abstract. A novel decomposition technique is used to extract the centrality dependence of 
di-jet properties and yields from azimuthal correlation functions obtained in Au+Au collisions 
at ^/s NN =200 GeV. The width of the near-side jet shows very little dependence on centrality. 
In contrast, the away-side jet indicates substantial broadening as well as hints for for a local 
minimum at A(f> = n for central and mid-central events. The yield of jet-pairs (per trigger 
particle) slowly increases with centrality for both the near- and away-side jets. These observed 
features are compatible with several recent theoretical predictions of possible modifications of 
di-jet fragmentation by a strongly interacting medium. Several new experimental approaches, 
including the study of flavor permutation and higher order multi-particle correlations, that 
might help to distinguish between different theoretical scenarios are discussed. 



1. Introduction 

QCD calculations on the lattice indicate a transition from hadronic matter to a deconfined phase 
of quarks and gluons (QGP) at extremely high temperature and energy density ^E]. Heavy ion 
collisions can deposit large amounts of energy in a small collision volume. Thus, they provide 
excellent pathways to energy densities and temperatures pertinent to such a phase transition. 
Indeed, several recent results from Au+Au collisions at the Relativistic Heavy Ion Collider 
(RHIC) indicate the creation of matter with energy density believed to be sufficiently high for a 
state of matter with non-hadronic degrees of freedom to be formed jHj. The dynamical evolution 
of this matter appears to be hydro dynamically driven, implying that the produced medium is 
strongly interacting [U 03 H3 |S1 E] . This matter has been observed to strongly suppress the 
number of high transverse momentum particles |1U1 1111 [T21 113j and to quench the away-side jet 
in central Au+Au collisions |14| . The suppression is commonly ascribed to radiative energy loss 
of hard scattered partons traversing the high energy density matter prior to the formation of 
hadrons [£3 HB1 El HH1 OH H3 • 

The influence of parton-medium interactions on jet properties such as the number of jet- 
associated partner particles per trigger hadron, and jet-topologies, remains an open question 
of great current interest. Consequently, much recent effort is targeted at understanding how 
jets couple to the strongly interacting medium produced in energetic Au+Au collision at RHIC 
j2U E2 E3 El E3 E2 E3 EE1 EH • At the center of this discussion is the possible creation of jet 
induced "conical flow" or "mach shocks" analogous to a sonic boom in a fluid [25 j l26l l2?l l28j. 



If observed, such a phenomena could provide the means for reliable estimates of the speed of 
sound in the nuclear collision medium produced at RHIC. 

In this contribution, we utilize azimuthal angular correlation functions for charged hadrons to 
examine the influence of medium-effects on jet characteristics and yields. The use of correlation 
functions provides a novel and effective probe for such effects because they provide simultaneous 
access to both the (di-) jet signals and harmonic contributions arising from hydrodynamic 
evolution of the system. 



2. Data Analysis 

The analysis presented here uses Au+Au data (y / SAw=200 GeV) provided by RHIC in the 
second running period (2001) and recorded by the PHENIX collaboration. The full PHENIX 
detector setup is described elsewhere |3()j . Charged tracks for this analysis were reconstructed 
in the two central arms of PHENIX, each of which covers 90 degrees in azimuth. Tracking was 
performed via a drift chamber and two layers of multi-wire proportional chambers with pad 
readout (PC1.PC3)|30|. A combinatorial Hough transform in the track bend plane was used 
for pattern recognition 31 . Most conversions, albedo and decays were rejected by requiring a 
confirmation hit within a 2 a matching window in the PC3. Collision centrality was determined 
using the beam-beam counters (BBC) and zero degree calorimeters (ZDC) by cutting in the 
space of BBC versus ZDC analog response 32 . 

Following the approach commonly exploited in HBT analyses', we define the azimuthal 
angular correlation function as the ratio of a foreground distribution N cor , constructed with 
correlated particle pairs from the same event, and a background distribution N m i x , obtained by 
randomly pairing particles from different events within the same multiplicity and vertex classes. 

Here, A<p = (</>i — 4>2) is the azimuthal angle difference between particle pairs formed with one 
hadron from a high-pr "trigger" bin (hereafter labeled A) and another hadron from a lower pt 
selection (hereafter labeled B). A particular asset of constructing a correlation function is that 
azimuthal acceptance and detector efficiencies automatically cancel in the ratio of foreground 
to background distribution. The area normalized correlation function then gives the probability 
distribution for detecting correlated particle pairs per event within the PHENIX pseudorapidity 
acceptance (\r]\ < 0.35). 

Figure ^ shows the centrality dependence of di-hadron correlation functions obtained from 
Au+Au collisions with px selections of 2.5 < px < 4.0 GeV/c and 1.0 < pt < 2.5 GeV/c for 
trigger and associated hadrons respectively. The correlation function for the most peripheral 
event sample (cf. Fig. U) shows two distinct peaks; a narrow peak at A<p = and a somewhat 
broader peak at A<p = tt. This asymmetic pattern is compatible with the very well known 
azimuthal patterns expected for di-hadron correlations resulting from di-jet fragmentation. In 
contrast to the peripheral events, the correlation functions for more central collisions show less 
prominent jet-driven asymmetries and appear to be dominated by large anisotropies which can be 
linked to harmonic correlations resulting from elliptic flow. It is noteworthy that the correlation 
functions for central - mid-central collisions show a minimum below that expected for harmonic 
contributions (A(j) = vr/2). Such a shift strongly suggests an away-side jet which is significantly 
broader than in peripheral collisions. 

A particularly important feature of the correlation functions shown in Figure ^ is the 
observation that they can be fully accounted for if one assumes that only two sources of 
correlations contribute to the measured azimuthal correlation function. That is, a (di-)jet 
correlations and and harmonic correlations due to elliptic flow j!41 1331 1341 135j . 



Figure 1. Correlation 
functions, C(A(f>), for trig- 
ger particles in 2.5 GeV/c< 
Pt <4.0 GeV/c and 
associated parti- 
cles in 1.0 GeV/c< 
<2.5 GeV/c. The 
bands show the harmonic 
contribution within the 
systematic uncertainty. 

3. Decomposition of jet and harmonic contributions 

It follows from the above discussion that a thorough investigation of possible modifications to the 
(di-) jet pair distribution in Au+Au collisions requires the separation of the jet-signal from the 
underlying harmonic "background" . Reliable procedures for such a decomposition are detailed 
in Refs. |361 I37| . Consequently, only the main steps of the procedure are outlined below. 

It can be shown j^H] that the pair correlations from the combination of flow and jet sources 
are given by 

C^(A0) = a o [C^(A0)] + J(A0), (2) 
where C^ B (Ac/)) is a harmonic function of effective amplitude V2, 

Cfi B (A0) = [1 + 2v 2 cos2(Act>)}; v 2 = (v^ x ). (3) 

and J(A(j)) is the (di-)jet function. By rearranging Eq. [21 one gets 

J{Acp) = C AB (A0) - a C£ B (A^). (4) 

J(A</>) can be extracted if the normalization (a D ) can be constrained and the harmonic amplitude 
of the background (V2) can be reliably determined. To fix the value of a D we assume that the (di- 
)jet function has zero yield at the minimum (ZYAM) A<j) m i n . This is equivalent to demanding 

min)- However, no explicit or implicit assumption for the functional 
form of J(A0) is required. For the results presented here, the minimum is determined by 
requiring the background to coincide with a functional fit to the correlation function. We obtain 
V2 by measuring the single particle distributions relative to the reaction plane t/^r- The reaction 
plane is reconstructed in the beam-beam counters and corrected for dispersion following the 
procedures outlined in Ref. [101 • An important aspect in this regard is the large (pseudo) rapidity 
separation (A \r]\ > 2.75) between each BBC and the PHENIX central arms. This relatively large 




rj separation is believed to substantially reduce possible jet contributions to the determination 
of V2 [HE] • 

Before applying the ZYAM decomposition technique to the analysis of data, we have tested 
its reliability via a detailed set of Monte Carlo simulations that accounted for the <f) and the rj 
acceptance of PHENIX (S3- 

After applying the ZYAM condition to the PHENIX data and utilizing the measured values 
of i)2 (t>2 = xv 2))i the harmonic contribution to the correlation function is given by the solid 
bands in Fig. ^ The jet-pair distribution is then obtained by building the difference between 
the correlation function and this harmonic contribution. It can be shown, that the integral of 
this distribution is proportional to the fraction of jet-correlated particles per event [HHIISZI- This 
fraction can be extracted by building the ratio of the sum of J(A</>) and the sum of C(A</>) (over 
all bins in A</>), 

EiC(A0i) w 

One readily obtains the efficiency corrected pairs per trigger or conditional yield CY from this 
fraction via multiplication by the ratio of the average number of detected particle pairs per 
event (N^), to the product of the detected singles rates (N^), (N^), followed by a final 
multiplication with the efficiency corrected sing les rate {Nf ff ), for the lower p T bin B P>1I37]. 

iN AB \ 

4. Results 

Figure 13 shows the centrality dependence of the extracted jet-pair distributions normalized to the 
number of trigger particles. The most peripheral bin exhibits an azimuthal pattern compatible 
with regular vacuum fragmentation of (di-)jets. That is, a near-side peak at A<p = and an away- 
side peak at Acf) = ir. The jet-pair distributions obtained at other centralities reveal additional 
striking features. In order to quantify these features, we separate the jet-pair distributions into 
two regions and assign Acf) = — A(j) m i n and A^ m j n — Acfr = tt to the near- and away-side jets 
respectively. Both pieces of the jet-pair distribution can then be characterized by their widths 
(RMS) and their respective per-trigger yields. The resulting values so obtained are plotted in 
Fig. |21 as a function of the number of participants. Similar results for the 0-20% most central 
d+Au collisions are also included in Fig. |3]to facilitate a comparison between the Au+Au and 
d+Au systems. 

Let us first follow the centrality evolution of the near-side jet. The shape of this jet is 
essentially unchanged across centralities (cf. Fig. [2J). However it's per-trigger yield rises with 
increasing collision centrality suggesting some degree of modification [by the medium] which 
does not have a strong influence on the near-side jet shape. In contrast to the near-side jet, the 
shape of the away-side jet-pair distribution broadens significantly as one moves from peripheral 
events to mid-central events (cf. Fig. In addition, the mid-central distributions (20-40% and 
10-20%) indicate an away-side jet distribution that exhibits a local maximum at Acp = 2ir/3 
and an apparent minimum at Acp = tt. These observations are consistent with a picture that 
involves little or no quenching of jets in peripheral events and strong centrality dependent 
medium modifications of the away-side jet shape in the mid-central and central events. It is 
interesting to note that both the near- and away-side per trigger yields (cf. Fig. |3J) show a mild 
rise with increasing centrality, suggesting a possible medium induced modification of both near- 
and away-side jets. The prominent "hump-backed" shape of the away-side jet is consistent with 
recent conjectures of a strong coupling between such jets and the high energy density matter 
that they traverse [25 ^ l2fil l2?| l2~4] . 



Figure 2. ZYAM- 
subtracted jet-pair distri- 
butions l/N trig dN/d(A4>). 
The dashed (solid) his- 
tograms show the distribu- 
tions that would result from 
subtracting the harmonic 
contribution after increas- 
ing(decreasing) the product 
V2 = V2 x V2 by one inter- 
val of the systematic error. 
The hatched area indicates 
the systematic error on the 
background determination 
(see text). 

It is instructive to study the influence of systematic variations in v% on the jet-pair 
distribution. The solid (dashed) lines in Fig|2 indicate the conditional yield distributions that 
is extracted after subtraction of a v<i product that is lowered (raised) by one interval of the 
systematic error respectively. The systematic error on v<i is dominated by the uncertainty on 
the reaction plane dispersion. The dotted line shows the jet-pair distribution resulting from 
a subtraction of a v% product lowered by twice the systematic uncertainty. In the latter case, 
the local minimum at A(f> = tt is no longer significant. However, the away-side peaks for all 
centralities smaller than 60%, remain significantly broader than in the most peripheral event 
sample. 

5. Summary of charged hadron results 

To summarize, we have used a novel decomposition technique to extract absolutely normalized 
jet pair distributions from correlation functions. Our results indicate jet pair distributions for 
peripheral Au+Au collisions which suggest very little, if any, medium induced modification of 
jet properties. By contrast, the away-side jet measured in central and mid-central collisions 
show strong modifications resulting from interactions with the medium. At this juncture, it is 
very interesting that the observed jet modifications are in qualitative agreement with several 
theoretical conjectures of how jet characteristics can be modified by a strongly interacting QCD 
liquid [21 [23 [23 123 [23 

6. Towards the Future 

Future detailed studies of the modification of jet characteristic by the strongly interacting QCD 
medium will require new measurements as well as enhancements in our analysis techniques. In 
what follows, we give a few examples which show that such developments are already underway. 
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Figure 3. (a) Per trigger yields for near- and away-side peaks in the jet pair distribution; and 
(b) Widths (RMS) of the peaks. The triangles denote results for 0-20% most central d+Au 
events from a recent PHENIX analysis. 



6.1. Extraction of the jet function via harmonic extinction 

For mid-central collisions, Fig. ^ shows pair correlations with strong harmonic components 
Cfj B (A(f>) and relatively weak jet components J(A0). Only after subtraction of C^ B (A0) is 
the true shape revealed for J(A^) (cf. FigEl). It would be desirable to select a data set that 
contained only those correlations due to the jets. In recent work we have developed a new 
technique which allows data selection such that the harmonic correlations are extinguished ((V2) 
= 0) [HZj- This technique utilizes A(j) distributions where the trigger particles are selected in a 
window perpendicular to the reaction plane |1J . For these distributions V2 can change phase and 
it is possible to select a cut angle 4> c so that V2 = 0. In Fig. |3]we demonstrate the extinguishing 
of harmonic correlations. The filled circles show the inclusive A(f> distribution obtained for a 
simulation with an input x>2 ~ 0.16. The filled squares show the out-of-plane A<f) distribution 
after the cut angle 4> c , is set to the extinction value, ie. <j) c = 4> x t |37j . The latter distribution is 
flat and demonstrates that this technique gives a good method to determine the jet correlations 
directly from a data set, without the blurring effect from harmonic correlations mediated by the 
reaction plane. The latter technique provides clear advantages for future jet studies at RHIC. 

6. 2. Flavor permutations 

The study of jet fragmentation and the influence of possible interactions of the jet with the 
medium in Au+Au collisions will undoubtedly require investigations which seek to map out the 
flavor dependence (ie. the dependence on particle species) of these effects. Correlation functions 
measured with fully identified particles at intermediate and high px should serve as a powerful 
tool to unravel different aspects of the interaction between jets and the medium. Recently, 
we have begun to measure such correlation functions at RHIC. Figure El shows preliminary 
flavor permutation correlation functions for the centrality range 20-40%. The trigger particle 
Pt bin spans 1.5 < pr < 2.5 GeV/c and the assoicated particle has been selected from the 
range 0.8 < pt < 1-5 GeV/c. One can observe strong anisotropies and asymmetries [in 
these correlation functions] reminiscent of those discussed above for charged hadrons. More 
importantly, one can clearly see that these anisotropies and asymmetries strongly depend on 
the flavor permutation (ie. the the identity of each particle which is used in the pair correlation 
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Figure 4. Simulated cor- 
relation function for uncon- 
strained particles (filled cir- 
cles) and for a trigger par- 
ticle constrained within the 
cut angle <p c = (p xt per- 
pendicular to the reaction 
plane (filled squares), see 
text. The results are for 
a pure harmonic simulation 
with t>2 ~ 0.16. 



measurement). For example, the baryon-baryon correlation function shown in Fig. El (a) does not 
appear to exhibit any asymmetries and shows a clear minimum at A(p = ir/2, ie. the correlation 
function is essentially harmonic. This suggests little, if any, contributions from (di-) jet sources 
to the di-baryon correlation functions at these transverse momenta. On the other hand, the 
meson-meson correlation function shown in Fig. [5] (c) is relatively asymmetric, suggesting that 
jet fragmentation products are predominantly mesons. These data suggest that a detailed study 
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Figure 5. Correlation 
functions for identified trig- 
ger particles (1.5 < px < 
2.5GeV/c) and identified 
partner particles (0.8 < 
Pt < 1.5GeV/c). The cen- 
trality selection is 20-40%. 



of correlation functions as a function of flavor permutation holds much promise for providing 
detailed insights into jet fragmentation and its modification in a strongly interacting medium. 



6.3. Higher order multiparticle correlations 

In recent work we have also been pursuing analyses' which involve higher order multi-particle 
correlations. Such correlations have been shown to carry additional topological information 
compared to a two particle correlation analysis |l2]. Additional shape information from higher 
order correlations could potentially provide a powerful tool for distinguishing between different 
scenarios for jet-modification by the nuclear medium. For example, "mach cone" like structures 
would potentially manifest themselves in a rather distinctive pattern when correlating more 
than two particles. Such analyses for identified and unidentified particle multiplets are currently 
being pursued with vigor. 
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